Calculation of the positron bound state with the copper atom 
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Abstract 



A new relativistic method for calculation of positron binding to atoms is 
Qh| presented. The method combines a configuration interaction treatment of 

the valence electron and the positron with a many-body perturbation theory 
description of their interaction with the atomic core. We apply this method 
jrt ' to positron binding by the copper atom and obtain the binding energy of 170 

c/3 , meV (±10%). To check the accuracy of the method we use a similar approach 

to calculate the negative copper ion. The calculated electron affinity is 1.218 
►|j^.^| eV, in good agreement with the experimental value of 1.236 eV. The problem 

^ ■ of convergence of positron-atom bound state calculations is investigated, and 

means to improve it are discussed. The relativistic character of the method 
and its satisfactory convergence make it a suitable tool for heavier atoms. 
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I. INTRODUCTION 

Bound states of positrons with neutral atoms have not been detected experimentally 
yet. For a long time the prevailing view was that neutral atoms do not bind positrons. 
For example, Aronson et al |[l|] proved that positron binding to hydrogen is not possible, 
and Gertler et al [0] showed that a ground-state helium atom could not bind a positron. 
In a number of calculations positron binding was observed for alkalis and second column 
atoms 0-§[. However, important physical effects, such as virtual or real positronium (Ps) 
formation, were neglected in those works. As a result, the binding was largely considered 
as an artifact of the approximations used, or the positron bound states found were unstable 
against Ps emission. This situation has clearly changed now. Firstly, a many-body theory 
calculation by Dzuba et al. |^ indicated that atoms with larger dipole polarizabilities and 
ionization potentials greater then 6.8 eV (Ps binding energy) can bind positrons, and pre- 
dicted positron binding energies for Mg, Zn, Cd and Hg. Subsequently, a number of recent 
calculations have shown and even proved, for a few lighter atoms, that positron-atom bound 
states do exist, [^0. 

For the problem of positron-atom binding the atoms should be divided into two groups: 
those with the ionization potential / smaller than 6.8 eV, and those with / > 6.8 eV. For the 
former the lowest fragmentation threshold of the positron-atom system is that of a positive 
ion and a Ps atom. Consequently, positron binding to such atoms should rather be described 
as binding of the Ps to the corresponding positive ion. Indeed, the 'ion + Ps' component 
in their wave function is large, as shown by the calculations for Li-e"'', Na-e"*" and He2'^S'-e"'" 
||7|j9Hl2|. For atoms with / > 6.8 eV the positron-atom bound state is indeed an 'atom + e"*"' 



system, at large positron-atom separations. However, the process of virtual Ps formation 
in this system is very important 0, especially when I is close to 6.8 eV. This effect makes 
positron-atom bound states a strongly correlated atomic system. The correlations in it are 
stronger than those one finds in its electron analogues, atomic negative ions. This feature 
makes the positron-atom bound complexes very interesting for the atomic theory. This also 
makes them a challenging testing ground for applications of modern numerical methods of 
atomic structure calculations. 

The main difficulty in calculations of positron interaction with atoms comes from the 
strong electron-positron Coulomb attraction which leads to virtual positronium formation 
0). One can say that it gives rise to a specific short-range attraction between the positron 
and the atom, in addition to the usual polarizational potential which acts between a neutral 
target and a charged projectile p^THl^. This attraction can not be treated accurately by 
perturbations and some all-order technique is needed. In our earlier works ||6| JT^JT9[] we 
used the Ps wave function explicitly to approximate the virtual Ps-formation contribution 
to the positron-atom interaction. The same physics may also explain the success of the 
stochastic variation method in positron-atom bound state calculations (see []I2| and Refs. 



therein). In this approach the wave function is expanded in terms of explicitly correlated 
Gaussian functions which include factors exp(— ar^?^) with interparticle distances Tj^. Using 
this method Ryzhikh and Mitroy obtained positron bound states for a whole range of atoms 
with both / < 6.8 eV (Li, Na, and He 2^5), and / > 6.8 eV (Be, Mg, Zn, Cu, and Ag). This 
method is well suited for few-particle systems. Its application to heavier systems is done by 
considering the Hamiltonian of the valence electrons and the positron in the model potential 



of the ionic core. However, for heavier atoms, e.g., Zn, the calculation becomes extremely 
time consuming [|l5l, and its convergence cannot be ensured. 

Another non-perturbative technique is configuration interaction (CI) method widely used 
in standard atomic calculations. This method has been applied to the positron-copper bound 



state in |T^. In this work the single-particle orbitals of the valence electron and positron are 
chosen as Slater-type orbitals, and their interaction with the Cu"*" core is approximated by 
the sum of the Hartree-Fock and model polarization potentials. The calculation shows slow 
convergence with respect to the number of spherical harmonics included in the CI expansion, 
Lmax = 10 being still not sufficient to extrapolate the results reliably to Lmax -^ oo. 

In the present work we calculate the ground states of Cu-e^ and Cu~ systems using a 
CI calculation within a spherical cavity of finite radius R. This procedure facilitates the 
convergence of the CI expansion in the difficult positron-atom case, and we show how to 
extrapolate the results to the R ^ oo limit. The CI method which we use is based on the 
combined relativistic configuration interaction and many-body perturbation theory method 



(CI+MBPT) developed in our earlier work pO| for precise calculations of many-electron 



atoms with more than one valence electron. It was shown there that correlations between the 
core and valence electrons are very important and often contribute more to the energy than 
the correlations between the valence electrons. The core-valence correlations are included 
into the effective CI Hamiltonian of valence electrons by means of many-body perturbation 
theory. This allows us to achieve high accuracy in calculations of atomic energies and 
transition amplitudes. In the present work we adapt this approach to the positron problem. 
As a single-particle basis for the CI calculations we use S-splined [^ Hartree-Fock wave 



functions in the cavity of finite radius R. The B-spline technique has been successfully used 



in atomic calculations for many years (see, e.g., review ||22|) and has been recently incor- 
porated with the CI+MBPT method ||2^. The use of 5-splines ensures good convergence 
of the CI calculation with respect to the number of radial orbitals. Convergence is further 
controlled by varying the cavity radius, while the effect of a finite cavity size on the energy 
of the system is taken into account analytically. 

We have chosen the copper atom for the positron bound-state calculations for several 
reasons. First, this atoms looks like a good candidate for positron- atom bounding. It has a 



large polarizability of 40 a.u. |2^, and its ionization potential / = 7.724 eV |2^ is not too 
far from the Ps binding energy of 6.8 eV, which ensures a sizable contribution of virtual Ps 
to the positron-atom attraction. Second, copper has a relatively simple electronic structure 
with only one valence electron above closed shells. This makes the positron-copper problem 
effectively a two-particle problem well suited for application of the CI+MBPT method. 
Third, there are accurate experimental data and a number calculations for the energy of the 
copper negative ion. Thus, we can test our method on Cu~ and compare the results with 
those obtained by other techniques. Last but not least, the existence of the positron-copper 
bound state was predicted by Ryzhik and Mitroy |0 in the framework of the stochastic 



variational method, which allows us to compare the results obtained with the two different 
techniques. 

II. METHOD OF CALCULATION 



A. Effective Hamiltonian 

We use the relativistic Hartree-Fock method in the V'^'^ approximation to obtain the 
single-particle basis sets of electron and positron orbitals and to construct an effective Hamil- 
tonian. The main point for this choice is the simplicity of the MBPT, as discussed in Ref. 
pi| . The self-consistent potential is determined for the Cu"*" ion and the single-particle 
states of the external valence electron and the positron are calculated in the field of the 
frozen core. 

The two-particle electron-positron wave function is given by the CI expansion, 

^(re, Vp) = Y, C,j^i{re)(f)j{rp), (1) 

id 

where ipi and (pj are the electron and positron orbitals respectively. The expansion coef- 
ficients Cij are to be determined by the diagonalization of the matrix of the effective CI 
Hamiltonian acting in the Hilbert space of the valence electron and the positron, 

-'^eflf ^ /^e + "-p + "-ep; 

h^ = cccp + (/3 - l)mc^ + Vr' - VZ:n' + ^e, 

hp = cap +{(3- l)mc' + V/"^ + Sp, (2) 

'^ep 1 f ^epi 

I e '-pi 

where he and hp the effective single-particle Hamiltonians of the electron and positron, 
and hf,p is the effective electron-positron two-body interaction. Apart from the relativistic 
Dirac operator, he and hp include the direct and exchange Hartree-Fock potentials of the core 
electrons, V^~^ and V^~j^, respectively. The additional S operators account for correlations 
involving core electrons (see [^ for a detailed discussion). We calculate S using the second- 



order MBPT in the residual Coulomb interaction. Eg describes the interaction between 
the valence electron and the electrons of the core. All four second-order diagrams for the 
Eg are presented in Fig. |I[ Sp is the correlation interaction between the positron and the 
core. In the second-order Ep is represented by a sole digram in Fig. ^j. Both operators are 
often called correlation potentials, because these non-local operators can be included into 
the equations for the single-particle orbitals together to the Hartree-Fock potential. Eg and 
Ep are energy-dependent operators, which are different for the electron and the positron. 
They are calculated separately for each partial wave, (si/2, Pi/2, P3/2, etc.). However, at 
large distances both operators have the same asymptotic behaviour, 

Ee(r,r'), Ep(r,r')^-|^5(r-r'), (3) 

where a is the dipole polarizability of the atomic core. This asymptotic form comes from the 
dipole contribution of the first diagram in Fig.^ for the electron, and diagram in Fig.^ for 
the positron. Formula (^ with some empirical cut-off at small distances is often used as an 
approximation for the correlation potentials, and is usually called 'polarization potential'. 



Sep is another type of correlations between the external particles and and core electrons. 
It can be described as screening of Coulomb interaction between the external electron and 
positron by the core electrons. There are in all six second-order diagrams for Egp. Three of 
them are shown in Fig. ^. The other three can be obtained from them by mirror reflection 
with respect to the vertical axis. When the electron and the positron are well outside the 
atomic core Egp is given by the following asymptotic expression, 

S.p(r..rp)=.^^^. (4) 

e p 

Similarly to Eq. (||), this formula is often used to construct rough approximations for Sep- 
Such potentials are called 'di-electronic correction', or 'two-body polarization potential'. 

Diagrammatic expansions in Figs. [^, H and |^ enable one to include valence-core corre- 
lations in an ah initio manner. To increase the accuracy of the calculations higher-order 
contributions to S can be taken into account effectively, by introducing a numerical factor 
before E. For example, the coefficient for Eg can be chosen by fitting the energies of the 
neutral atom states to the experimental data. In doing so the important non-local structure 
of the operators is preserved. 



B. Basis set 



We use 5-spline basis functions |21[ to calculate the diagrams for S and to construct 
the single-particle orbitals for the CI expansion ([l|). For this purpose the atomic system is 
confined to a cavity of radius i?, and the wave functions are set to zero aX r = R. For a suffi- 
ciently large R the error introduced by this boundary condition is very small for atomic-size 
binding energies, ~ exp(— 2/T;i?), where k is related to the binding energy as e^ = K^fi^ /2m. 
However, for weakly bound states, e.g. those of the positron with the atom, this error has 
to be considered more carefully (see below). The interval [0, -R] is divided into a number 
of segments and 5-splines are constructed on them as piecewise polynomials of a certain 
degree. They are bell-shaped overlapping smooth functions. With an appropriate choice of 
the radial mesh they can approximate atomic wave functions to a very high precision. Note 
that it is not convenient to use S-splines directly in CI or MBPT calculations because of 
their non-orthogonality. Instead, we use their linear combinations which are eigenstates of 
the single-particle Hartree-Fock Hamiltonian. This ensures orthogonality, allows to separate 
core and valence states and improves convergence, since only a relatively small number of 
lower Hartree-Fock eigenstates are sufficient for the convergence of the CI calculation. This 
also means that while we use the same 5-splines for the electron and positron states the 
resulting single-particle basis states are different, because the Hartree-Fock Hamiltonians 
for the electrons and positrons are different. Another advantage of the use of 5-splines is 
that the convergence can be controlled by the cavity radius R (its reduction leads to a more 
rapid convergence), while its effect on the energy is taken into account analytically. 



C. Effect of finite cavity size 

The choice of the cavity radius R (see above) is dictated by a compromise between the 
convergence rate and the required accuracy of the calculations. On one hand, the radius must 
be large enough to accommodate the wave function of the state under investigation, e.g., 
the positron-atom bound state. On the other hand, smaller radii mean faster convergence, 
both with respect to the number of radial orbitals and, which is especially important for 
positron-atom calculations, to the number of angular harmonics. This effect is very strong 
since convergence is determined by the cavity volume which is proportional to i?^, and having 
a smaller radius means that one needs fewer basis states to describe the wave function. 

The problem of convergence is crucial for the positron-atom interaction. As discussed 
in the Introduction, the positron tends to form virtual Ps with the external atomic electron 
)| JT7|JI8[] . The positronium radius rpg ~ 2ao can be small compared to the characteristic 



size of the positron-atom bound state wave function, r ~ 1/ k ^ Oq, where oq is the Bohr 
radius. To describe Ps at large separations from the atom expansion (|T]) needs to be ex- 
tended to very high values of angular momentum L and principal quantum number n to 
account accurately for the virtual Ps formation. This problem is well known in positron- 
atom scattering calculations, see e.g. ||2^. Smaller cavity radii force virtual Ps to be at 



smaller distances, thereby improving the convergence significantly. However the energy of 
the system is affected. Therefore, the convergence and the accuracy of the calculation can 
be really improved only if the effect of a finite-radius cavity on the energy is taken into 
account. 

To consider the effect of cavity on the energy of the system let us consider the problem 
of a particle weakly bound in an s state by a finite-range potential. 'Weakly bound' here 
means that the binding energy is much smaller than the typical scale of the potential. This 
is definitely true for positron-atom bound states whose binding energy is much smaller than 
1 eV. To determine the radial wave function x(r) at large distances it is sufficient to impose 
on it a boundary condition 



Xdr 



(5) 



at the outer radius r = a of the potential well ||2^. The k parameter is related to the 
energy of the bound state e = —k^H /2m, and determines the asymptotic form of the wave 
function, x(r) ~ Ae~'^^ . 

The boundary condition is unchanged when we place the system in the cavity of finite 
radius R, R > a, provided the energy of the bound state is still small. However, the wave 
function must now turn into zero at the cavity radius, x{R) = 0. This shifts the energy 
of the weakly bound state up from e to some other value e/j, which depends on the radius 
of the cavity. The Schrodinger equation for a < r < R, where the potential is vanishingly 
small, is 

|;|^ + ^«M-o, (6) 

After solving it with boundary conditions (|^) and x{R) = 0? o^is obtains a negative eigen- 
value, €r = —nj^fi^ /2m, where 

6 



K = kr/ tanh[/T;ij(i? — a)], (7) 

if R is not too small, R — a > k,^^. As one can see, for R —^ oo the solution of Eq. (^, 
kr, approaches its asymptotic value k, and the energy in the cavity e/j — > e. For a smaller 
cavity radius the eigenvalue becomes positive, Er = kj^h'^ /2m, where kji is found from 

K = k^/ taii[kR{R — a)]. (8) 

This means that the state which is bound may appear as unbound due to the effect of the 
cavity. Equation (H) is valid for kji{R — a) < |. Otherwise, k < 0, and the energy is too 
high, so that it remains positive even when the cavity wall is removed. 

Equations (|^ and (H) can be used to find the infinite-cavity energy e = —H?h /2m from 
the energy 6^ calculated for the finite cavity radius R. It is important that these formulae 
are insensitive to the detailed shape of the atomic potential, and depend only on the atomic 
radius a. The value of a can be estimated from the position of the classical turning point 
Tc, in the potential for an external atomic electron, 

Tc 2aoz/2 ' 

where v is the effective quantum number of the valence electron. Beyond the turning point 
the valence electron's wave function decreases exponentially, as exp(— r/aoz/). Therefore, a 
resonable estimate for is 

a = rc + aou = (2z/^ + z/)ao (9) 

For copper (/ = 0.28349 a.u., z/ = 1.33) this gives a ~ Soq. A more accurate value of a can 
be found by applying Eq. (|^) to two bound-state calculations performed with two different 
cavity radii R. The uncertainly in the value of a is in fact unimportant, as long as we 
consider weakly bound states for which na <^1. 

Note that the wave function is also affected by the finite cavity size. This should be 
taken into account in calculations of the annihilation rate and other matrix elements. The 
annihilation rate is proportional to the probability of finding the positron close to the atom. 
For a -C -R the wave function at r ^ a is affected via normalization only. The change of the 
normalization can be found by comparing the normalization integral for r > a calculated 
numerically within the cavity, 

rR 
J a 

with the analytical value 

/■oo 1 



a 



2k 



III. RESULTS AND DISCUSSION 

A. Copper negative ion 

To test the method and find out what accuracy can be achieved we first apply it to the 
copper negative ion. This is an effective two-particle problem technically very similar to 
the positron-copper interaction considered above. It should be mentioned that for Cu~ only 
the electron S-operator is involved (Fig. |I|), and for the screening of the electron-electron 
interaction, instead of the diagrams on Fig. ^, one must use similar diagrams presented in 
20| (Fig. 4). The results of calculations for Cu and Cu~ are presented in table | together with 



the experimental values. The energies are given with respect to the Cu"*" core. The accuracy 
of the Hartree-Fock approximation is very poor. The binding energy of the 4s electron in 
neutral Cu is underestimated by about 20%, while the negative ion Cu~ appears altogether 
unbound (its energy lies above that of the neutral atom). The inclusion of core- valence 
correlations (S) does improve the energy of the neutral atom, but the negative ion is still 
not bound. The standard CI method, in contrast, takes into account the valence- valence 
correlations, while neglecting the core-valence correlations. It does produce binding for the 
negative ion, but the binding energy is almost two times smaller than the experiment value. 
Only when both core-valence and valence-valence correlations are included the accuracy 
improves significantly. It is equal to 2.6% for the ionization potential of the neutral atom and 
10% for the electron affinity, which is quite good for a relatively simple ab initio calculation. 
The remaining discrepancy is mostly due to third and higher-order correlation corrections 
in E, since the configuration expansion for Cu~ converges rapidly, and the corresponding 
error is small. 

To simulate the effects of higher-order terms in S and thus further improve the accuracy 
of calculations we introduce numerical factors before the Eg operators to fit the lowest s, 
p and d energy levels of the neutral copper atom. These factors are fs = 1.17, fp = 1.42 
and /d = 1.8 in the s, p and d channels, respectively. Table |l| shows that these factors also 
significantly improve the calculated electron affinity. It is natural to assume that the same 
procedure should work equally well for the positron-atom problem. 

Results of other calculations of the electron affinity of copper are presented in table 0. 
Note that only a coupled-cluster (CC) method produces a result more accurate than ours. 
It is interesting to mention among other results the results by Mitroy and Ryzhikh p!3|JT6[ 



who calculated Cu~ for the same purpose as we do, i.e., to gauge the accuracy of their 
method for the positron-atom problem. Their first result for electron affinity to copper, 
0.921 eV, was obtained by the stochastic variational method, while another result 0.916 eV 
was achieved in the standard CI calculation. Both methods are variational in nature and 
differ basically by the form of the trial two-electron wave function. Since the two results 
agree well with each other, good convergence has probably been achieved in both methods. 
However there is a noticeable discrepancy between their result and the experimental electron 
affinity value. From our point of view the most probable source of this discrepancy is 
approximate treatment of the S operator of the valence-core interaction. In their works 
Mitroy and Ryzhikh use approximate expressions for the core polarization potentials, based 
on asymptotic formulae (|^) and (^), which include only dipole core polarization in the local 
energy-independent form. Note again that the actual S operator is energy dependent. It is 



different for different angular momenta, and for the electron and positron, while approximate 
expression (^ is always the same. Note also that the screening operator Sep depends not only 
on the states involved but also on the multipolarity of the Coulomb integral. Approximate 
formula (|) describes the dipole part of screening only, however, other Coulomb multipoles 
are also screened. Even though the largest contribution to screening comes from the dipole 
term, monopole and quadrupole screening can not be neglected. For example, monopole 
screening directly contributes to the diagonal Hamiltonian matrix elements in important 
configurations like 4s^ in Cu~, while dipole screening affects only the off-diagonal matrix 
elements. 



B. Positron binding to copper 

The binding energy of Cu~ is about 0.045 a.u. It corresponds to a bound-state parameter 
K, ^ 0.3, and the cavity does not have a noticeable effect on the calculated energies of Cu~, 
let alone Cu. The relative error introduced by the cavity can be easily estimated from Eq. 
(0), and even for a moderate R = 15ao it does not exceed 0.1% for the electron affinity. 

For the positron bound state the situation is different. As indicated by the calculation of 
Ref. |]T3[, the k value for the Cu-e"*" bound state is about 0.1. This is why we have performed 



the calculation of the positron-atom bound state using two different cavity radii, R = ISoq 
and R = 30ao, to make sure that convergence is really achieved. The convergence pattern 
with respect to the number of basis states used is illustrated in Figs. ^ and ||. Both plots 
show the energy of the electron-positron pair moving in the field of Cu"*", with respect to the 
energy of the Cu"*" ion (in atomic units). Empty circles correspond to R = ISag, while solid 
ones correspond to R = SOag. Dashed line shows the ground-state energy of the neutral 
copper atom. The positron-atom state is bound when its energy is below the dashed line. 
Fig. ^ shows the electron-positron energy of Cu-e"*" as a function of the number of radial 
basis functions in each electron and positron partial wave, n. The total number of partial 
waves is fixed by Lmax = 10. Note that convergence is visibly faster for the smaller cavity 
radius. For R = ISoq saturation begins at n ^= 10 while for R = 30ao the same level of 
saturation can be seen only at n ~= 18. Fig. ^ shows the Cu-e''' energy as a function of 
the number of partial waves included, while the number of radial wave functions in each 
wave is fixed at ra = 16 for R = 15ao and n = 22 ioi R = SOoq. Saturation can be clearly 
achieved for both radii at Lmax ~ 10. The difference in energy at the last (lowest) points 
for R = 15ao and R = 30ao in both figures is the effect of a finite cavity radius. It shifts 
the energy obtained in the R = 15ao calculation up with respect to the R = 30ao result. 
This effect can be easily taken into account using the formulae presented in section [11 C| . It 



turns out that the results for both cavity radii coincide, i.e., yield identical k from Eq. (^, 
for the atomic potential radius of a = 5.5ao- The final binding energy obtained is 0.0062 
a.u., or 170 meV. This should be compared to the result of Ryzhikh and Mitroy [|1^], which 
is 0.005518 a.u. or 150 meV. From the discussion of the accuracy of calculations which 
follows we conclude that the difference between two results is within the accuracy of both 
methods. A similar value is achieved in the CI calculation [0, which used 14 radial orbitals 



in each partial wave up to Lmax = 10, after extrapolation to Lmax = oo. However, the latter 
procedure has considerable uncertainties. 

There are several factors which affect the accuracy of our calculations. 



The accuracy of calculation of S and contributions of higher-order correlations. This 
can be estimated by comparing calculations with and without the fitting parameters, 



as discussed in section |IIl7\| . The introduction of the fitting parameters for the electron 
part of the correlation operator Eg reduces the binding energy by about 0.0009 a.u. 
However, the relevant uncertainty must be considerably smaller. Firstly, we saw that 
the use of fitting parameters really improves the calculated electron affinity of copper. 
We should expect the same effect for the positron binding energy. Secondly, the effects 
of the fitting parameters on the electron and positron operators Eg and Sp largely 
cancel each other. 

• Incompleteness of the basis set. We have seen from Figs ^ and ^ that the level of 
convergence achieved is very high and the corresponding uncertainty is small. Never- 
theless, there is a hidden uncertainty related to the radial coordinate mesh used, the 
number of splines and other parameters which determine the details of the numerical 
procedure. Varying these parameters shows that their effect on the binding energy 
does not exceed 0.001 a.u., when estimated conservatively. 

• Finite cavity radius. This effect on the binding energy calculated ai R = SOa^ is very 
small (~ 0.0001 a.u.). Since the results for R = 15ao and R = 30ao coincide for very 
reasonable value of the positron-atom potential radius a = 5.5ao, it is reasonable to 
believe that the corresponding uncertainty is very small too. 

Note that the difference between our calculated electron affinity of copper and the exper- 
imental value is 0.00066 a.u. If this value is compared with the numbers presented above, it 
is evident that it also gives a reasonable estimate of the accuracy of the calculation of the 
positron-copper binding energy (about 10%). Since the accuracy of calculations by Mitroy 
and Ryzhikh [J13l is not discussed in their paper the only thing we can do to estimate it is 



to use the same approach. Their best result for Cu~ differs from the experimental value by 
0.0116 a.u. (20% of the positron binding energy). If we adopt this value as the uncertainty 
of their result for the positron binding by copper, we see that the two results for Cu-e"'" 
bound state agree with each other within the accuracy of the methods. 
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TABLES 

TABLE I. Ground state energies of Cu and Cu~ calculated using different approximations (a.u.). 

Cu 

RHF^ -0.23830 

RHF + i:^ -0.27672 

Cr -0.23830 

CI +T,'^ -0.27672 

CI + / X S^ -0.28394 

Experiment^ -0.28394 

'^Relativistic Hartree-Fock; a single-configuration approximation, no core-valence correlations are 
included. 

Single-configuration approximation, core-valence correlations are included by means of MBPT. 
■^Standard CI method. 

CI+MBPT method, both core-valence and valence-valence correlations are included. 
^T, for s-wave is taken with factor / = 1.18 to fit the Cu ground state energy. 
^References [^,33]. 



Cu- 


Electron affinity 


-0.20309 


-0.03521 


-0.27280 


-0.00392 


-0.26424 


0.02594 


-0.31802 


0.04130 


-0.32869 


0.04475 


-0.32935 


0.04541 



TABLE II. Electron affinities of Cu (eV). Comparison with other calculations and experiment. 

Affinity Ref. Method 

Theory 
1.06 [||] Nonrelativistic MR CI calculations 

1.01 [p3| MR CI calculations in the DK no-pair formalism 

1.199 g 

1.236 Ql Relativistic coupled cluster method 

0.921 [^ Nonrelativistic stochastic variational method 

0.916 H] Nonrelativistic CI method 

1.218 Present work 

Experiment 
1.226 
1.2358 
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12 3 4 

FIG. 1. Second-order diagrams for the self-energy of the valence electron (Sg operator). Sum- 
mation of excited electron states a and j3 and core hole states m and n is assumed. 



a P b 




a 



FIG. 2. Second-order diagram for the positron self-energy (Sp operator). Double line denotes 
positron states. 
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FIG. 3. Screening of the positron-electron Coulomb interaction (Sgp operator) 
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FIG. 4. Energy of Cue"*" as a function of the number of radial electron and positron basis 
functions in each partial wave (Lmax = 10). Open circles are for R = 15ao, and solid ones for 
R = 30ao. 
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FIG. 5. Energy of Cue"*" as a function of maximal orbital momentum of the electron and 
positron orbitals in the CI expansion. Open circles are for R = 15ao, and solid ones for R = 30aQ. 
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